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In this paper, we consider fermionic systems in discrete spacetime evolving with a strict notion 
of causality, meaning they evolve unitarily and with a bounded propagation speed. First, we show 
that the evolution of these systems has a natural decomposition into a product of local unitaries, 
which also holds if we include bosons. Next, we show that causal evolution of fermions in discrete 
spacetime can also be viewed as the causal evolution of a lattice of qubits, meaning these systems 
can be viewed as quantum cellular automata. Following this, we discuss some examples of causal 
fermionic models in discrete spacetime that become interesting physical systems in the continuum 
limit: Dirac fermions in one and three spatial dimensions, Dirac fields and briefly the Thirring 
model. Finally, we show that the dynamics of causal fermions in discrete spacetime can be efficiently 
simulated on a quantum computer. 



I. INTRODUCTION 

The idea that there is a maximum speed of propagation 
of information has become one of the fundamental prin- 
ciples of physics. In particular, it appears in relativistic 
quantum field theories, the most interesting example of 
which is the standard model, which provides a unified 
framework for describing the effects of the strong, weak 
and electromagnetic forces. In this paper we are going 
to look at discrete spacetime quantum systems with the 
requirement of strict causality. Quantum systems in dis- 
crete space but continuous time with local Hamiltonians 
do not generally satisfy this notion of causality. (For an 
example illustrating why local Hamiltonians in discrete 
space typically lead to instantaneous propagation of in- 
formation, see appendix [Aj) There are a few reasons that 
causal quantum systems in discrete spacetime are inter- 
esting. 

The first is simulation: simulations of physics usually 
start by discretizing continuous degrees of freedom. Of 
particular interest to us is the simulation of relativistic 
quantum field theories by a quantum computer. Simu- 
lation of quantum physics is likely to be the first prac- 
tical application where a quantum computer could out- 
preform a classical one [T] . The idea that quantum com- 
puters may be better equipped than classical computers 
to simulate quantum systems dates back at least as far 
as Feynman in [2] . The basic idea was that quantum sys- 
tems evolving in a local way ought to be efficiently sim- 
ulable on a quantum computer using local operations, 
something that appears impossible on a classical com- 
puter. This was shown to be true for continuous time 
quantum systems with local Hamiltonians in [T] and, for 
fermions in particular, in [3]. These approaches rely upon 
breaking up the total Hamiltonian H into a sum of k 
body terms Hi, where k does not grow with the system 
size, and using a Trotter expansion in terms of e - lff i'/™ 
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to recover e _ *. For causal unitaries, however, it is not 
obvious how to simulate the evolution efficiently. We 
show how to do this for fermions by decomposing the to- 
tal unitary evolution operator exactly into local unitaries 
on qubits, such that the evolution can be efficiently sim- 
ulated on a quantum computer. 

The second reason we may be interested in discrete 
spacetime quantum models that are strictly causal is that 
nature itself may be discrete, and it is plausible that 
even at the smallest length scales there is a strict form of 
causality. The idea that the spacetime continuum breaks 
down in some way at very small length scales is often 
embraced with a view to constructing a theory reconcil- 
ing gravity with quantum mechanics, such as causal sets 
[I]. What is really interesting is that there are causal 
discrete spacetime models that become interesting rel- 
ativistic quantum field theories in the continuum limit. 
We describe some of these in section [iVj It is encouraging 
that, although these discrete systems do not have Lorentz 
symmetry, it is recovered in the continuum limit. Cur- 
rently, the only such models are fermionic; a particularly 
interesting example is given in [5 , which becomes the 
Thirring model in two dimensional spacetime. Ideally, 
one would like a general recipe giving a discrete space- 
time model that converges to a given quantum field the- 
ory in the continuum limit if it exists. Taking such limits, 
however, is a complicated process as couplings must be 
renormalized. 

The final reason causal quantum models in discrete 
spacetime are interesting is that they provide discretized 
models of relativistic systems that are well defined. De- 
spite its remarkable successes in explaining physics at 
very high energy scales, quantum field theory in the con- 
tinuum has yet to be put on a firm mathematical footing. 
Here, because we work with a discrete lattice, we have 
a regulator, so the infinities that plague quantum field 
theory do not appear. 

The models we will look at are quantum systems in 
discrete spacetime with a strict notion of causality. They 
evolve over each timestep via a causal unitary opera- 
tor, which essentially means that localized observables 
can only spread a finite amount over one timestep. Our 
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main focus will be fermionic systems: one of our re- 
sults is the proof of a general principle (theorem [I] in 



section |III A ) that says that causal fermionic evolution 
can be decomposed exactly into a product of local com- 
muting fermionic unitaries, which is analogous to a re- 
sult of [6] for quantum cellular automata. This extends 
their maxim that "unitarity plus causality implies local- 
izability" to systems of fermions. We also extend this to 
systems of fermions and bosons that may be interacting. 

Causal fermionic systems in discrete spacetime are 
fermionic analogues of Quantum Cellular Automata 
(QCA) [BHS]. These are discrete spatial lattices that have 
finite quantum systems (with associated finite dimen- 
sional Hilbert spaces) at each spatial point. QCA evolve 
in discrete time via a causal unitary operator [25] . A use- 
ful picture to have in mind is a spin lattice that evolves 
over discrete timesteps in a causal way. QCA are inter- 
esting for many reasons. They are universal for quantum 
computation |9] , with the particularly nice property that 
they are implcmcntablc by applying local unitary gates 
O [7] . We show that the dynamics of causal fermionic 
systems in discrete spacetime can be viewed as subsec- 
tors of the dynamics of QCA (theorem [2b . This is anal- 
ogous to a result in jTDl E] , which maps local fermionic 
Hamiltonians to local spin Hamiltonians. Their results 
for local Hamiltonians provide a way of thinking about 
fermions without the need for anticommuting operators 
on different sites. Theorem [2] extends this to systems 
of causal fermions in discrete spacetime. It is interesting 
that in both cases anticommuting fermionic operators are 
not necessary to describe the dynamics. 

The breakdown of this paper is as follows. We start 
in section [TT] by going through the properties of fermions 
and the Jordan- Wigner Transformation, which allows us 
to represent fermionic operators by operators on qubits 
(two dimensional quantum systems, which we can think 
of as spin 1/2 particles). In section III we look at causal 
fermionic evolution in discrete spacetime and show that 
it always has a decomposition in terms of local fermionic 
unitaries. We also point out that this extends to possi- 
bly interacting bosonic and fermionic modes. In section 
|III B| we show that we can view any causal fermionic 
evolution as a causal discrete time evolution of qubits, or 
equivalently as a subsector of the evolution of a quantum 
cellular automaton. Next, in section [TV] we give some dis- 
crete spacetime models that become interesting systems 
in the continuum limit: first, we reproduce the evolution 
given in [5j [12l [13] of discrete spacetime Dirac fermions 
in one dimensional space and its corresponding local im- 
plementation on qubits. Second, we look at a similar 
discrete fermionic model that behave like Dirac fermions 
in three dimensional space in the continuum limit, orig- 
inally given for a single particle in [2]. In section IV D 



we look at the representation of the Dirac field in discrete 
spacetime. This is followed by a discussion in section [V] 
of the implications of these results for simulation. After 
concluding remarks, we extend these results to infinite 
lattices via C*-algebras in the appendix. 



Note that we set h = c = 1 throughout. 



II. FERMIONS AND THE JORDAN- WIGNER 
TRANSFORMATION 

For simplicity, here we will restrict the set of spatial 
points to be a finite-sized d-dimensional lattice, so that 
positions are labelled by vectors of integers. To introduce 
translational symmetry, it will sometimes be helpful to 
treat this lattice as a torus (i.e. to introduce periodic 
boundary conditions). The results all have natural ex- 
tensions to different geometries, but it will be useful to 
have this particular example in mind for when we dis- 
cuss locality and causality. We postpone a discussion of 
systems on infinite lattices to appendix [F] as these in- 
volve additional complications that are not necessary to 
understand the main ideas. 

Now suppose that we have fermionic modes at each 
point. A natural example of this is a system where each 
site can be occupied by spin up or spin down electrons. 
We denote the state with all modes unoccupied by |£1). 
Then we define creation and annihilation operators at^ 
and ag^, where x labels the position and \i labels the 
extra degree of freedom at each site. These satisfy the 
canonical anticommutation relations: 



(11.1) 



where Sgg = 1 if x = y and is zero otherwise. We also 
have that a^lfi) = 0, and the Hilbert space is spanned 
by all possible products of creation operators ett^ acting 
on For example, the state at^atjO) has a fermion 
at x and a fermion at y, with internal degrees of freedom 
/i and v respectively. 

An extra requirement that we make of systems of 
fermions is that physical observables are not only self- 
adjoint but are linear combinations of products of even 
numbers of creation and annihilation operators. In par- 
ticular, in continuous time systems, all physical Hamil- 
tonians have this property. For example, the Hubbard 
Hamiltonian is 



H = —a V~Yat a,-,, + at a-,, 



<xy> M 



(II.2) 



where [i G {t, 1} labels the extra degree of freedom (spin 
in this case), (xy) denotes nearest neighbour pairs, and 
a, U > are real parameters. The first term describes 
electrons hopping, while the second term is an on site 
Coulomb repulsion. 

Because fermion creation and annihilation operators 
anticommute regardless of the spatial separation between 
them, these operators are, in a sense, non local objects. 
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If we want to represent a fermionic system by a system of 
qubits, we will have to take this into account. For exam- 
ple, suppose we want to represent a line of N fermionic 
modes with no internal degree of freedom and positions 
labelled by x £ {0, 1, N — 1} by TV qubits. It is natu- 
ral to take |00...0) to represent the state with no fermions 
present. Now we can choose the representation of a J on 
the qubits to be 



Al = \(X -iY ) = |l)o{ 



(II.3) 



where X, Y and Z are the standard Pauli operators 
[25] . The subscript implies that these operators act on 
all other qubits like the identity, meaning, for example, 
X n = X (g> / ® I...® I. Because a\ and a\ anticommute, 
we cannot simply pick a\ to be represented by A[. In- 
stead, we can satisfy the anticommutation relations by 
taking 



,t = 



4nx 



(II.4) 



y<x 



This is known as the Jordan- Wigner Transformation |15j . 
It preserves the anticommutation relations because of the 
strings of Zs. In fact, we are free to choose the order- 
ing in the product above however we want. To see this, 
we give a more general Jordan- Wigner Transformation, 
which is particularly useful for higher dimensional lat- 
tices. Given N fermionic modes, we associate a qubit to 
each mode. Next, we assign a unique label to each site, 
tt(x) £ {0, N — 1}, and define 



4=4 

7r(y)<7r(z) 



(II.5) 



which also satisfies the anticommutation relations. For 
the special case of a line of fermions with n(x) = x we re- 
It is sometimes useful to use different 
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cover equation 
ordering schemes for different problems. 

If there are multiple fermionic modes at a lattice site 
(for example, due to spin, or different types of parti- 
cle) , then the qubit representation will include a separate 
qubit at that lattice site for each mode. Then the order- 
ing tv(x, /i), where [i distinguishes different modes at site 
S, will assign a unique number to each mode. It is natu- 
ral to choose the ordering such that ir(x, /J,) for the set of 
modes at each site are consecutive. This means that even 
products of fermionic creation and annihilation operators 
at the same site are local in the qubit representation. 
With the ordering in (II.4I, local fermionic operators 



on a line are represented by local operators on the cor- 
responding line of qubits. In higher spatial dimensions, 
however, it is not generally true that local fermionic op- 
erators correspond to local qubit operators. In fact, even 
for a ring of fermions this is not generally true. Luck- 
ily, there is a way of getting around this that involves 
introducing auxiliary fermions, given in [101 lllj and re- 
produced in appendix [D] which we will use to extend the 
results of the following sections to systems of fermions in 
arbitrary spatial dimensions. 



III. CAUSAL FERMIONS 

In the following sections, where we prove our main re- 
sults, we look at discrete time systems. In continuous 
time, it is natural to work with a Hamiltonian as it de- 
termines the evolution via the Schrodinger equation. In 
a discrete time picture, however, there is no Schrodinger 
equation, so it is more natural to work directly with the 
unitary operator that acts on the state each timestep. 

To define notions of locality and causality, it is help- 
ful to define the neighbourhood of a spatial point. Here 
we take this to mean the set of points that are at most 
1 lattice step away in each spatial direction (taking into 
account the periodic boundary conditions if necessary) 
|27) . The neighbourhood of a point is therefore a d- 
dimensional hypercube of side length 3, centered on that 
point. Note that this definition of the neighbourhood is 
not critical to the proofs: for example, if there were next 
nearest neighbour interactions, we could consider larger 
hypercubes. 

We say that a fermionic operator is localized on a spa- 
tial region R if it can be written in terms of creation 
and annihilation operators corresponding only to R. Be- 
cause of this, and the fact that they must be sums of even 
products of creation and annihilation operators, localized 
observables from non overlapping regions of space always 
commute. 

Next, we define a causal fermionic unitary. Note that 
we work in the Heisenberg picture. 

Definition 1. A fermionic unitary U is causal if, for 
any x and fi, Wag^U is localized in the neighbourhood of 



The definition of a causal unitary ensures that over one 
timestep information cannot propagate farther than one 
step in each spatial direction. Note that we can construct 
any operator localized on a region from creation and an- 
nihilation operators corresponding to that region. So, in 
particular, this definition ensures that local observables 
do not spread very far after one timestep. 

In continuous time systems, the Hamiltonian of a sys- 
tem of fermions is a sum of even products of creation 
and annihilation operators. In particular, this implies 



that 



-iHt 



commutes with the annihilation operator b if 



H does not contain any terms with b or b> . We take it for 
granted that the discrete time dynamics of fermions we 
consider also have this property. This means that, given a 
system of fermions evolving via U, we can add additional 
fermionic modes whose creation operators anticommute 
with the original fermion creation and annihilation oper- 
ators while commuting with U. (We could just assume 
that U = e tA , where A is a sum of even products of cre- 
ation and annihilation operators, but this is less general 
and not useful for infinite systems.) 

Finally, we have the following useful lemma, which is 
proved in appendix [B| 

Lemma 1. The inverse of a causal fermionic unitary is 
also a causal fermionic unitary. 



4 



A. Local Decomposition 

Here we will give a local unitary decomposition for 
causal fermionic unitaries analogous to that given for 
QCA in [3 [7]. In short, we will show that the full unitary 
evolution can be decomposed into a product of unitaries, 
each of which is localised in a particular neighbourhood. 

Constructing this local decomposition requires us to 
consider the joint evolution of the system of fermions 
together with an identical copy of that system. Let us 
denote the creation operators for the original system of 
fermions by at and the creation operators for the cor- 
responding modes of a copy system by ot . 

It will be useful to define a fermionic swap operator. A 
unitary implementing the fermionic swap at O bg , 
meaning Sg^ag^Sg^ — bg^ and Sg^bg^Sg^ = ag^, is 

= exp[i»(6^ - aV){b 3li - a Sv )}. (HI.l) 

See appendix [C] for more details about this operator. 

The local decomposition of a causal evolution of 
fermions is derived in the following theorem. 

Theorem 1. Given a finite system of fermions with cre- 
ation operators at evolving via the causal unitary Ua, 
the evolution of two copies of this system via UaU b , 
where U B is equivalent to Ua but acting on the o|, 
fermions, can be decomposed into local fermionic uni- 
taries: 

u A u B = n( s *>) n^^^], (in.2) 

where U B Sg u U B are commuting local fermionic uni- 
taries. 

Proof. First, 

II( S ^) UPbS^uI] = SUbSUI, (III.3) 

where 

S = TJ(S^) (III.4) 

is the unitary that swaps all modes. And, because S 
swaps all modes, St/gS = Ua- It follows that 

II( S *>) II^S^l] = U A Ul (III.5) 

x,[i y,v 

Furthermore, U B Sg^U B is a local fermionic unitary be- 
cause Sg^ is just 

ex P[*|(4 M ~ 4p)( 6 *M ~ °*/*)]> ( IIL6 ) 

and so 

U B S S M = exp]i^(b% - 4 M )(% - a^)], (III.7) 



where b' s = U B bg^ U B , which must be localized within 
the neighbourhood of x because U B is causal. Hence 
U B Sgfj,U B is also localised within the neighbourhood of 
x. (Naturally, we are thinking of each mode labelled by 
x as being at the same site as its copy, which is why the 
unitaries UbS$ij.U b are local.) The fact that the different 

unitaries J/bS^C/Jj commute follows from [Sj M ,S^] 
' □ 

This theorem tells us that the joint causal evolution of 
two copies of a system of fermions can be decomposed 
into a product of local unitaries. Now note that, given 
any state of the physical system of fermions and its 
copy, it follows that, for any measurement operator M.a 
on the physical fermions, 

(1>\UbU\MaUaUIWi = (1>\U^MaUa\1>), (III.8) 

so this joint evolution of both the physical and auxil- 
iary fermions is just as good as the original evolution 
but with the advantage of being decomposable into local 
fermionic unitaries. For example, \ip) could be any state 
of the physical fermions with all auxiliary modes unoc- 
cupied, such as ;^§(<4i + a y- u )\tyi where is the state 
annihilated by all ag^ and bg^. 

In some cases there is a natural local unitary decom- 
position of the evolution without the need to include a 
copy of the system: see the example in section [IV A| In 
general, however, this is not true. A simple counterex- 
ample is the unitary that shifts everything one step to 
the right every timestep. 

As an aside, note that a similar decomposition exists 
for systems of bosons and fermions, which may be inter- 
acting. To make the notation simple, we suppose that 
there is only one fermionic and one bosonic mode at each 
point, though the result still holds with multiple modes 
at each site. We denote by §g the bosonic swap operator 
between the mode at x with creation operator ct and its 
copy with creation operator dt. 

S s = cxp[i|(4 - cl)(dg - eg)}. (III.9) 

The local decomposition is encapsulated in the follow- 
ing corollary. 

Corollary 1. Given a finite system of fermions and 
bosons evolving via the causal unitary U A , the evolution 
of two copies of this system via U A U B , where U B is equiv- 
alent to Ua but acting on the copy system, can be decom- 
posed into local unitaries: 

U A UI 

= no 3 *) ij( s *-) n^ s ^] uius^uii (iilio) 

w x y z 

where U B SgU B and U B E>gU B are commuting local uni- 
taries. 
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B. Representation by Qubits 

Now that we have a local decomposition of a causal 
fermionic unitary, we can look at its representation in the 
qubit picture. In this and the following section, we will 
assume there is only one fermionic mode per site to make 
the notation simpler. (The extension to more than one 
mode per site is straightforward.) We assign qubits to 
the fermionic mode at x and its copy (created by at and 

bt respectively) . For now we will work with a finite line of 
points, so that with the natural choice of ordering for the 
Jordan- Wigner Transform (equation II. 4) the operators 
Sg are local unitaries in the qubit representation. Next, 
recall that UsSgUg is just 



exp[z-(&g-4)(&^ 
tfrrt 



(III.ll) 



-iHs 



Now note that b 1 = UsbgU B must be a linear combina- 
tion of odd powers of creation and annihilation operators 
on the neighbourhood of x. (This is proved in lemma [2] 
in appendix B ) So Ub&xU b has the form 



where 

Hg is a self adjoint operator localized on the neighbour- 
hood of x containing only even products of creation and 
annihilation operators. But with the natural ordering 
for the Jordan- Wigner Transformation in equation |II.4| 
Hg is localised on the neighbourhood of x in the qubit 
representation also. This means that in the qubit repre- 
sentation UbSsU b is localised on the neighbourhood of 
x. 

Hence, we can view this causal evolution of fermions as 
a causal evolution of qubits. In higher spatial dimensions 
or lines with periodic boundary conditions, we need to do 
more to ensure that our causal fermionic evolution can 
be represented by local unitaries acting on qubits. Wc 
see how to do this in the next section. (Note that the 
main points about simulation and models of quantum 
field theories can be understood without going through 
the details of this.) This proves the following theorem. 

Theorem 2. Any causal fermionic evolution in discrete 
spacetime is equivalent to a subsector of the causal evo- 
lution of a system of qubits, which is a type of quantum 
cellular automaton. 



C. More than one spatial dimension 

Now that we are considering higher spatial dimensions, 
in general there is no choice of ordering for the Jordan- 
Wigner Transformation such that all local unitaries in 
the fermion picture are local in the qubit picture. Nat- 
urally, we choose the ordering so that fermionic modes 
and their copies, which are associated to the same site, 
are consecutive in the ordering scheme. This means that 
each fermionic swap Sg term in the decomposition is still 
local in the qubit representation. 

The U B SgU B terms in the decomposition are more 
problematic. They are not necessarily localized unitaries 



in the qubit picture. Fortunately, we can circumvent this 
problem but at the price of adding auxiliary fermions, 
which also means that we need more qubits to represent 
this larger system of fermions. 

First, we already know that UbSsU b is a local 
fermionic unitary. And, as we saw in the previous sec- 
tion, it has the form e~ lHs , where Hg is a self adjoint 
operator localized on the neighbourhood of x containing 
only even products of creation and annihilation opera- 
tors. Next, we use an idea from jTDl [IT]. (For brevity, 
here we will just give a rough idea of how this works; see 
section|D]for full details.) Suppose, for example, that Hg 
contains the term agbg, where y is in the neighbourhood 
of x. Suppose also that this term is non local in the qubit 
representation: this will be because of strings of Z op- 
erators arising from the Jordan- Wigner Transformation. 
Now we introduce two new fermionic modes, one at x and 
one at y, with annihilation operators eg and eg. Because 
we choose the Jordan- Wigner Transformation ordering 
such that modes at the same site are consecutive, ageg 
and Cjjbjj are local in the qubit picture. This is because all 
Z operators in the qubit representation of ag correspond- 
ing to sites other than x are cancelled by those from the 
qubit representation of eg. 

Next, we make the replacement 



agb$ -t ag(im s mg)bg, 



where 



eg + c\ 



m V ~ C V + C y> 



(111.12) 



(111.13) 



which we can think of as Majorana fermions. The new 
term on the right hand side of 1 2| is local in the qubit 
picture, and, acting on a +1 eigenstate of imgmg, 



ag(imgmg)bff 



^x^y- 



(III. 14) 



Similarly, for any other terms in Hg that are non local in 
the qubit representation, we can add additional fermions, 
so that, acting on +1 eigenstates of all the additional 
terms like imgmg that we add, the resulting operator is 
equal to Hg. 

This means that UbS$U b is equivalent to a local 
fermionic unitary on a larger system that is local in the 
qubit picture. So we have extended theorem [2] to any 
dimension. 

Note that the fermionic unitaries UbSsU b commute. 
After introducing additional fermions, the new unitaries 
Vg implementing UsSgU B on the qubits do not neces- 
sarily commute when the neighbourhoods on which they 
are localized overlap. Notice, however, that the order in 
which they are applied does not matter when acting on 
a +1 eigenstate of the Majorana pairs inigmg. Further- 
more, we can apply many Vg simultaneously, provided 
they act on non overlapping regions. For example, for 
a line, we can apply all V x with x mod 3 = first, fol- 
lowed by all V x with x mod 3=1, followed by all V x with 
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x mod 3 = 2. So for a line we need only three steps to 
implement every V x . 

On one hand, it seems that we could have just ex- 
panded the causal fermionic unitary U as a sum of prod- 
ucts of creation and annihilation operators and applied 
the trick of adding Majorana fermions directly to this. It 
is not clear, however, that we can do this and preserve 
unitarity. Applying the prescription to a Hamiltonian 
worked because we could always preserve self-adjointness. 
Instead, our approach was to derive an exact decompo- 
sition of causal fermionic unitaries into a product of lo- 
cal unitaries, each of which can be written in the form 
e~ lHfS , where each Hg has support only on the neighbour- 
hood of x. We then applied the trick of adding Majorana 
fermions to each Hg. 

It is important from the point of view of simulation 
and for the extension to infinite lattices that the number 
of additional fermions we need to add per site does not 
depend on M , the number of original fermionic modes. 

IV. CONSTRUCTING MODELS 

In the following sections, we take a constructive ap- 
proach and examine specific examples of causal discrete 
spacetime models that become interesting continuum 
models as we take a continuum limit. We also repre- 
sent the evolution of such discrete spacetime systems by 
products of local unitaries on qubits. Note that the re- 
alization of this for the one dimensional Dirac equation 
appeared in O E3 US] . 

A. Discrete Dirac Fermions in One Dimension 

In this section, we will look at fermions in discrete 
spacetime that obey the one dimensional Dirac equation 
in the continuum limit. We will use theorem [T] to decom- 
pose the evolution into a simple product of local fermionic 
unitaries, and then we will give a representation of this 
in terms of qubits. Both of these steps allow us to re- 
produce the evolution given in (T2l [T3l [16]. (A single 
particle in discrete spacetime that obeys the Dirac equa- 
tion in the continuum limit appeared in [17] . The details 
of this continuum limit were studied further in [TM2"U] .) 

For simplicity, suppose that our discrete space is fi- 
nite, with sites labelled by n 6 {0, ...,N — 1} with pe- 
riodic boundary conditions. And suppose there are two 
fermionic modes at each site, labelled by I and r. As 
we anticipate getting Dirac fermions in the continuum 
limit, let us denote creation operators by ijp n a , where 
n e {0, ...,N - 1} and a G {I, r}. Also define 

*■ = (&) ■ (IV1) 

We define the matrices /3 and a>i in this basis to be 

"-(!i).«-(J-°i)- 



It will be convenient to work in momentum space: 
p = 2itk/N is the discrete momentum, where k € 
{ — ^2^, ^2^}, and we take N to be odd. The mo- 
mentum creation operators are 

^ = ^E elP "<a- (IV.3) 
* n 

Let us also suppose that over each timestep the sys- 
tem evolves via the unitary U = WT . In the continuum 
limit, the unitary W will contribute the mass term in 
the Hamiltonian, and the unitary T will contribute the 
momentum term. First, T is a conditional shift that has 
the effect 

^ - - (t 1 ;,;) ■ <™> 

We can write T in terms of the discrete momentum op- 
erator. The operator that translates VVs.a one step to the 
right is 

exp(-£P a ), (IV.5) 

where 

p 

So we have 

T = exp(-i[P r - Pi]) = exp(-i ^ PV£«nM- (IV.7) 

p 

We define W to be 

W = exp(-iM^2 V'p^p)- (IV.8) 
p 

Now, to take a continuum limit of this, we embed the 
N spatial points into a line of length L (with periodic 
boundary conditions), with lattice spacing e = L/N. We 
must also let the number of timesteps r grow as e —¥ 0, so 
let t = re, with t constant. To get a sensible continuum 
limit, we set M = me, where m is a constant. Defining 
p = p/e and ^p p = ip pi and using Trotter's formula |21) . 
we get 

lim U t/e = exp[-i V ipUpa! + mP)ip p t], (IV.9) 

P 

where the sum is now over p — 2i:k/ L, with k G Z. So 
the continuum limit corresponds to particles evolving via 
the Dirac Hamiltonian in one spatial dimension: 

H D =J2 ^4(P Q i + m ^P- (IV-10) 
p 

Furthermore, we could take L to infinity to recover Dirac 
fermions on an infinite line. 

Now that we have found the continuum limit, let us 
return to the discrete time evolution. The conditional 
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shift part of the evolution (the unitary T in equation 
IV. 4 1 shifts ip ni i to the left and ipn^ to the right. But we 
can think of this as one system of ip n ^ fermions evolving 
via a shift to the left and a copy of that system, the VVi.r 
fermions, evolving via the inverse unitary: a shift to the 
right. This allows us to apply theorem [T] to see that this 
is equivalent to the fermionic swaps 



4>n,l "H> 4> n -l,r, 

at each n, followed by 

V-Vr <-» 1pn,l, 



(IV.ll) 



(IV.12) 



at each n. Applying the local unitaries in ( IV.11[ ) followed 
by the local unitaries in (IV.12 1 reproduces the original 
conditional shift in (IV.4). 



Note also that the part of the evolution that models 
mass is also a product of local unitaries, since the unitary 
W in position space is 



C. Discrete Dirac Fermions in Three Dimensions 

In this section, we construct a causal discrete space- 
time model that, in the continuum limit, becomes a sys- 
tem of fermions obeying the Dirac equation in three spa- 
tial dimensions. To do this, we first construct discrete 
fermions that obey the Weyl equation in the continuum 
limit, and then extend this to fermions obeying the Dirac 
equation in the continuum limit. Note that this evolution 
for a single particle was given in |14j . 

Suppose that our discrete space is finite with periodic 
boundary conditions and sites labelled by three compo- 
nent vectors n with components in {0, N — 1}. Sup- 
pose also that there are two fermionic modes at each site, 

labelled by a = 1,2. Let us define the corresponding 

t 



fermion creation operators by ^l a . Also define 



*Pn,2 



(IV.16) 



expHM tfnfrM = Hexpi-iM^lP^). (IV.13) We take to be the Pauli matrices: 



So the evolution operator U is a product of local uni- 
taries. 

In the next section, we represent this discrete fermionic 
system on qubits. 



B. Representation by Qubits 

We associate a qubit to each mode, such that the ba- 
sis states \l)ni and |0)„; correspond to the presence and 
absence of a left handed fermion at the point n, and simi- 
larly the states |l) nr and |Q) nr correspond to the presence 
and absence of a right handed fermion at the point n. 

We then represent the fermion creation operators us- 
ing the Jordan- Wigner Transformation, with the ordering 
7r(n, I) = 2n and n(n, r) = 2n + 1, so that 



n Zk - 

■K(k,i;)<TT(n,}j,) 



(IV.14) 



With this choice of ordering, each local fermionic unitary 
we found in the last section is local in the qubit picture 
with the exception of the swap ip ; O V'jv-i r across the 
periodic boundary. For example, the swapping operator 
in (IV.ll I for n ^ becomes 



exp[i-[A [ 



(n—l)r 



(IV.15) 



which is a local unitary on the qubits. The non-local 
swap crossing the periodic boundary can be dealt with 
be adding an extra pair of fermionic modes at positions 
and N — 1 and using the trick described in section |III C| 
Finally, as W is a product of on-site unitaries (equation 
IV.13 ), these will be local in the qubit picture. This form 



of the evolution of discrete Dirac fermions on a (infinite) 
line was given in [T3J, [T31 US] . 



1 

1 



CT2 = 



-i 

1 



^3 = 



1 

-1 



(IV.17) 



Again, it is simpler to work in momentum space: p = 
2nk/N is the discrete momentum, where the components 
of k take values in {— ^f^}, and we take N to 
be odd. The momentum creation operators are 




< n,a ' 



with 



(IV.18) 



(IV.19) 



via 



Suppose that over each timestep the system evolves via 
the unitary U = T1T2T3, where Ti are conditional shifts 
in each spatial direction: 



Ti = exp(— i^Pi^goW'p)- 



(IV.20) 



Note that each Ti is causal. Analogously to equation |IV.9| 
in section |IV A[ by applying the Trotter formula, these 
fermions obey the Weyl equation in the continuum limit. 
In other words, in the continuum limit they evolve via 
the Hamiltonian 



(IV.21) 



with the sum ranging over all p = 2?tk/L, where k has 
integer components. 



As in the one dimensional case (section IV A I, we can 
view the discrete evolution operator U = TiT 2 T 3 as a 
product of local unitaries. This is because each Ti is a 
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conditional shift operator in the ith spatial direction that 
can be rewritten as a product of local swap operations. 
We elaborate on this in appendix [E] 

Now, to get fermions obeying the Dirac equation in the 
continuum limit, we need four fermionic modes at each 
site. Call the creation operators for these modes xbt 

and ipt l , with a = 1, 2. In the continuum limit r and 
I will correspond to right handed and left handed modes 
respectively. We define 



"06 



p.r 



(IV.22) 



where each of the components of this vector has two com- 
ponents: 



*» " a,ld " (& 



^p>,2 



(IV.23) 



Let the evolution operator be U = WT, where in the 
continuum limit W will contribute the mass term, and T 
will contribute the momentum term in the Hamiltonian. 
Let T = T1T2T3, where Tj are conditional shifts acting 
differently on the I and r fermions: 



T t = exp(+i^2p i ip^ l (7 i ip$j) 
p 

X exp(-2 ^ Pii>p, r (T ii>P,r)- 
P 

Note that both terms commute. Define 

where / is the 2x2 identity matrix. Also define 



(IV.24) 



(IV.25) 



<Ti 
0~i ) ■ 



Then Tj can be rewritten as 

Ti = cxp(-i ^2 Pi^p^p)- 



(IV.26) 



(IV.27) 



And, similarly to the one dimensional case, we have 
the mass term 



expHM^^/%). 



(IV.28) 



As in the one dimensional case, in the continuum limit 
we get fermions evolving via the three dimensional Dirac 
Hamiltonian: 



(IV.29) 



with the sum ranging over all p — 2itk/L, where k has 
integer components. 



The unitary determining the evolution of these discrete 
Dirac fermions, U — WT1T2T3, is equivalent to a product 
of local unitaries. This is because W is a product of 
on-site unitaries mixing between ipft tT and tp^ i and each 
Tj can be decomposed separately into local unitaries for 
both r and I modes (as we saw for the Weyl case above) . 
To view this as a causal evolution of qubits, however, 
additional fermionic modes (and hence additional qubits) 
would have to be introduced to simulate the effects of 
fermionic anticommutation. Nevertheless, we saw how 
to do this in section ITlI CI 



D. Fermionic Fields 

Now let us turn to Quantum Field Theory (QFT), 
where in the usual approach fields are fundamental, and 
particles emerge after quantization. Take the Dirac field, 
with field operators ip a (x), where in three dimensional 
space a € {1, 2, 3, 4}, but in one or two dimensional space 
a G {1, 2}. The field operators obey the Schrodinger pic- 
ture anticommutation relations 1221 



{i/j a (x),ipp(y)} = 
{Mx)^l(y)} = 6 a p6 (3 \a 



y)- 



(IV.30) 



Ideally, we would like to represent individual electrons 
and positrons at site x by fermionic modes at x, but it 
is not clear how to do this. To see why, note that the 
continuum field operator in three dimensional space is 




(a^a s a (f)e^ + b s Jv s a (p)e-^% 



P 8=1,2 



(IV.31) 



where E p = +y / \p\ 2 + to 2 , a s J and bp create electrons 
and positrons with momentum p and spin state labelled 
by s € {1,2}, and u s a {p) and v s a {p) are four component 
eigenvectors of the Dirac Hamiltonian, satisfying 

a 

^2v s J(p)v^{p) = 2E p S sr (IV.32) 

a 
a 

See f22] for more details. Now the field operator in 



(IV.31 ) obeys the anticommutation relations (IV.30 1, but 



the electron part on its own, given by 
d 3 p 1 



A a (x) = 



a s pK(P)^ S , (IV.33) 



J P s=l,2 

does not satisfy 

{A a (x) : Al(y)} = 8 a0 6^(x-y). 



(IV.34) 
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This means we cannot assign separate spatial fermionic 
modes to represent both an electron field and positron 
field. This is related to the problem of not being able 
to localize electrons with only positive energy wavefunc- 
tions. So we will work with ip a {x) directly. 

In the continuum, the field operator ip a (x) obeys the 
Dirac equation. For a discrete model, we represent the 
field by fermionic modes at each point evolving via the 
discrete spacetime evolution that we considered in section 
IIV Al or IIV CI We know that the continuum limit of this 
is free fermions obeying the Dirac equation. This is es- 
sentially what is presented in (T3]. Each particle evolves 
in accordance with the Dirac equation in the continuum 
limit, but they do not always correspond to particles with 
positive energy. This is essentially because the vacuum 
is taken to be the state annihilated by "0s. a for all n and 
a. 

The naive vacuum defined by ip a (x)\ft) = for all 
x and a does not correspond to the physical vacuum 
in QFT, as the physical vacuum has no electrons or 
positrons present. Hence, we require that the physical 

vacuum is annihilated by all bt and at. We can cre- 

«/ p p 

ate this physical vacuum \Qd) by acting on |0) with all 
bp operators, which ensures that = 0. (This is 

equivalent to the Dirac sea picture. Viewed in that way, 
bp creates a negative energy electron, and bp creates a 
hole in the sea of negative energy particles.) So in the 
discrete case we need to consider a new vacuum state 
analogous to the Dirac sea state in the continuum. We 
postpone a detailed discussion of this problem for future 
work. 

Finally, note that in [S] a causal fermionic model in 
discrete spacetime is given that becomes the massive 
Thirring model in one spatial dimension in the contin- 
uum limit. In the massive Thirring model, evolution is 
governed by the Hamiltonian 



1 



H = H D + -g Uxf{x)jJx) 



(IV.35) 



where g is a constant, Hp is the free continuum Dirac 
Hamiltonia n in one spatial dimension (equation IV. 10 
in section 

i$ (x)a l ip(x 



IV A) and j^{x) is the current: j°(x 
ijj'(x)ip(x) and j l (x) 



V. SIMULATING CAUSAL FERMIONS 

We have seen how to represent a causal fermionic uni- 
tary by applying local unitaries to a lattice of qubits, 
which is a type of quantum cellular automaton. This 
tells us how to simulate the evolution of causal fermions 
on a quantum computer. 

From a complexity point of view, simulating the evo- 
lution of these systems can be done efficiently because, 
for Af fermionic modes, we have to apply 0(Af) local uni- 
taries: first the 0(Af) unitaries V£ M (the unitaries imple- 
menting Ub^>s^U\ in the qubit representation), followed 



by the qubit realisation of the 0(Af) fermionic swap op- 
erators. As we saw, we may need to include additional 
qubits to ensure that these operators are local in the 
qubit picture. Note also that the cost of applying Vgn 
does not grow with Af [25] . Furthermore, a lot of these 
operations can be done in parallel: all swap operations 
can be done simultaneously, while we can do many Vg^ 
operations at the same time, provided the areas on which 
these unitaries are localized do not overlap. As each Vg^ 
is localized on hypercubes with length of side 3 (because 
the evolution is causal), the time needed to implement 
one step of the evolution is 0(3 d ), where d is the lattice 
dimension, so the time does not depend on Af . 

So much for implementing the evolution on a quan- 
tum computer. We still have to prepare a +1 eigenstate 
of the pairs of Majorana fermions on the quantum com- 
puter. This can be done efficiently by using a method 
presented in [3] to deal with the strings of Z operators 
that arise in the qubit representation. The method is 
given in appendix [Gj 

Finally, there additional subtleties involved when try- 
ing to simulate quantum field theories. These are dis- 
cussed in |23j , which shows how to simulate 4> A theory on 
a quantum computer in a very different manner. 



VI. DISCUSSION 

Throughout this paper, we have discussed causal quan- 
tum systems in discrete spacetime. The first results we 
obtained were theorem [T] and corollary [T] which allow 
us to decompose causal unitaries into a product of local 
unitaries in a manner analogous to that of [6j[7]. Later, 
we used this and a method for mapping local fermionic 
Hamiltonians to local qubit Hamiltonians to prove theo- 
rem [2] which showed that causal fermionic quantum sys- 
tems in discrete spacetime can be viewed as lattices of 
qubits evolving causally in discrete time, meaning they 
are types of quantum cellular automaton. After dis- 
cussing specific discrete spacetime fermionic models, we 
showed why these systems can be efficiently simulated on 
a quantum computer. 

The next objective is to devise causal discretized mod- 
els that become interesting quantum field theories in the 
continuum limit, such as quantum electrodynamics or 
quantum chromodynamics, with the ultimate goal be- 
ing the construction of causal models that reproduce the 
entire standard model in the continuum limit. One of the 
reasons this is interesting is that this would allow sim- 
ulation of these models by a quantum computer, some- 
thing that seems feasible for purely fermionic systems. 
For bosonic systems any such program would require a 
cut-off to allow the state of each bosonic mode to be rep- 
resented by a finite number of qubits. 

We hope that the results contained here may help in 
some way to find causal discrete spacetime models that 
converge to interesting physical systems in the contin- 
uum limit. This would not only be useful for simulation 
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but also as mathematically sensible (discretized) quan- 
tum field theories with a strict notion of causality that 
offer an alternative view to current discretized models. 
Furthermore, the study of causal discrete models may 
even hint at physics beyond the standard model, particu- 
larly as it is sometimes suggested that at some small scale 
the notion of continuous spacetime may break down. 
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Appendix A: Continuous Time Models in Discrete 
Space 

Here we will see why we cannot construct interesting 
discrete space systems in continuous time that are strictly 
causal. Take a particle on a finite discrete line evolving 
in continuous time with some time independent Hamilto- 
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nian H. Our strict notion of causality would imply that, 
if the particle is at position at t = 0, then there is some 
T such that for t < T the particle has zero probability of 
being found outside a finite region R containing 0. 

This tells us that, given any position n outside R, we 
must have that, for all t < T, (n|e _lfft |0) = 0. Expand- 
mg e" 



-iHt 



(n\e- lHt \0) = (n| - iHt + O(t 2 )\0) = 
(n\iH\0) - (n\O(t)\0) = 0, 



(A.l) 



but the second term can be made arbitrarily small by 
taking t to be small, so the first term must be zero. Sim- 
ilarly, by looking at higher order terms in the expansion 
of e~ lHt , we see that (n\H l \0) = for any I. But this 
implies that (n\e- tHt \0) = for any t. It follows from 
this that the particle will not be found outside R at any 
t. 

Therefore, if the particle is ever going to propagate to 
a point ri, it happens instantaneously, though possibly 
with a very small amplitude. 



Appendix B: Proofs of Lemmas \T\ and [2] 

Here we prove lemmas [l] and [2] It will be useful to re- 
peat the extra requirement we made of causal unitaries in 
section [ITT] Given a system of fermions evolving via U, we 
can add additional fermionic modes whose creation oper- 
ators anticommute with the original fermion creation and 
annihilation operators while commuting with U . This is 
essentially a discrete time analogue of the requirement 
that Hamiltonians are sums of even products of creation 
and annihilation operators in continuous time systems. 

For simplicity of notation, we will assume here that 
these fermions have no extra degrees of freedom, though 
the extension to systems of fermions with extra degrees 
of freedom is straighforward. 

Lemma 2. Given a unitary U acting on fermions with 
annihilation operators ag, U'a^U is a linear combina- 
tion of odd products of fermion creation and annihilation 
operators. 

Proof. We write Wa s U = A odd + A even , where A odd are 
all the terms that are products of an odd number of cre- 
ation and annihilation operators and A even are all terms 
that are products of an even number of creation and an- 
nihilation operators. 

Our extra requirement above implies that we can add 
a fermionic mode with creation operator w , which anti- 
commutcs with all of the original creation and annihila- 
tion operators while commuting with U . But this implies 
that a s }U = {o f , A odd + A even } = 0, which is only 

possible if A even — 0. □ 

Next we prove Lemma [T] from section [TTT| 

Lemma [JJ The inverse of a causal fermionic unitary U 
is also a causal fermionic unitary. 



Proof. Lemma [2] tells us that WagU must be a linear 
combination of odd products of creation and annihilation 
operators. So, since U is causal, 



{U^a^U,a g } =0, 
and {U^aUl,a g } = 



(B.l) 



for all y when x is not in the neighbourhood of y. Then 

(B.2) 



{a ff , Ua^} = 0, 
and {at, Ua s U^} = 



for all y when x is not in the neighbourhood of y. But 
x is in the neighbourhood of y implies y is in the neigh- 
bourhood of x. It follows that Ua^W is localized on the 
neighbourhood of x, so W is causal. □ 



Appendix C: The Fermionic swap operator 

Given two fermionic modes with annihilation operators 
a and b, it is useful to define a unitary operator which 
swaps them: 

S^aS = b, S^bS = a. (C.l) 
Here we show that this unitary is given by 

(C.2) 



5 = cxp[i|(fe t -a f )(6-a)] 



which is self-adjoint. To see this, define two new 
fermionic modes 

c=^=(b-a), d=^=(b + a), (C.3) 



V2' 



V2' 



which satisfy the usual anticommutation relations. Then 
S — exp[i7rc^c], and 



SaS = e 1 ^ c ^- (d - c) e 
V2 



Wc = ^(d + c) = fe, 



where we have used the fact that (c^c)c = 0, and c(c^c) = 
c. Similarly, it is easy to see that SbS = a. 



Appendix D: Representing Local Fermionic 
Hamiltonians by Local Qubit Hamiltonians 

Here we will see the main idea from [TUJ [TT] . In section 



III C and appendix F 4 we apply this technique to local 
fermionic operators. In |10[ lllj. however, they apply it 
to local fermionic Hamiltonians, which are sums of local 
fermionic operators. 

To illustrate the idea, let us look at a local fermionic 
Hamiltonian. By including additional fermionic modes, 
we will construct a local fermionic Hamiltonian that 
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replicates the dynamics of the original Hamiltonian but 
is local in the qubit representation. As an example, take 



H = (4 a ?+ o }°«f)' 



(D.l) 



<xy> 



where (xy) denotes nearest neighbour pairs, and we 
are considering a rectangular lattice as shown in figure 
[I] With the Jordan-Wi gner Transform given by (|II.5[), 




FIG. 1: A lattice of 5 x 4 different sites, with Jordan- Wigner 
Transformation ordering, n(x), as shown. 

and the ordering for the Jordan- Wigner Transformation 
shown in figure [T] all of the vertical hopping terms in 
(D.l) will be non local in the spin picture. Furthermore, 



with any other choice of ordering, some terms in (D.l) 



would be non local in the spin picture. And, as we con- 
sider bigger and bigger lattices, the length of the strings 
of Z operators in the non local terms will grow. 

We want to eliminate these non local strings of Zs. To 
do this, we will introduce extra fermions to cancel the 
strings of Zs in the qubit representation. 

In particular, we introduce a pair of additional 
fermionic modes whenever two sites are connected in the 
Hamiltonian by a hopping term ata^ + ata^ that is not 
local in the qubit picture (for example, when tt(x) = 0, 
7r(y) = 5 in hgure [I]). We introduce one additional 
fermionic mode at site x and one at site y, with cre- 
ation operators denoted by cj- ^ and cj- s ^ respectively. 
The first index in the subscript gives the location of the 
ancillary fcrmion, and the second indicates the hopping 
destination. We define the operators 



m (x,y) = c {x,y) + C (g,g) 
m (y,x) = c {y,x) + c \y,g) ' 

which are self-adjoint and satisfy 

X^ix^y) : TTL(w,z) } = ^xw^yz 

{m{S,y) , a-z} = 0. 



(D.2) 



(D.3) 



These can be thought of as Majorana fermion operators. 



Now define the operator Mrg^gj 



lm (x,y) m (y,x) 



,7 in and 



note that it is self-adjoint and has eigenvalues +1 and 
-1 since Mf s g) = 1 and M {s ^ ^ -1. 
Next, we make the substitution 



atag + alas 



a x M (x,y) a v' 



(D.4) 



in the Hamiltonian in D.l whenever sites x and y are 
not adjacent in the ordering scheme. As all of the op- 
erators Mig^ commute, there exists a joint eigenstate 
with eigenvalue +1 for each operator. (This is the reason 
we must distinguish between the ancillary modes 
and C(^?). If we were to replace both with a single mode 

then the operators Mtg^ and M^ g^ would not com- 
mute.) When the ancillary modes are prepared in this 
state, then the action of the transformed Hamiltonian 
on the original fermions will be the same as that of the 
original Hamiltonian. 

It is natural to choose the new ordering such that 
fermions at the same site are consecutive because, as we 
mentioned in section[Hj this means that the product of an 
even number of fermionic creation and annihilation op- 
erators at the same site will always be local in the qubit 
representation. Then, because m^g^ is a fermionic oper- 
ator at site x, the operator atm^^ is local in the qubit 
representation. Therefore it follows that 



a l M (x,y) a y + a l M (x,y)ax 



(D.5) 



is local in the qubit representation. The same trick 

will work for any local fermionic Hamiltonian that is 

quadratic in fermion creation and annihilation operators. 

Furthermore, even if the fermionic Hamiltonian has 

higher order terms (provided they only contain even 

products of creation and annihilation operators), like 

a \ a \ a w a z + ata^affag, for example, we can still make 

this local in the qubit picture by adding more Majorana 

fermions. For this example, we can replace atata^ag by 

x y 



a x M (x,y) a^a-wM^^ ag. 



(D.6) 



By using this trick, any local fermionic Hamiltonian 
has a corresponding local qubit Hamiltonian. Also, each 
term in the qubit Hamiltonian connects the same sets 
of lattice sites as the corresponding term in the original 
fermionic Hamiltonian 

Furthermore, because [at, M($a] = for any x, y and 
z, it follows that we can act on a +1 eigenstate of all the 
^(S,y) terms with physical fermion creation operators at 
and still have a +1 eigenstate. Therefore the original 
fermionic dynamics can be viewed as a subsector of the 
dynamics of the corresponding local qubit Hamiltonian. 

Note that the requirement that the Hamiltonian is a 
sum of even products of creation and annihilation oper- 
ators was necessary to use this trick. Furthermore, the 
fact that we could preserve the self-adjointness of the 
Hamiltonian was also crucial. 

Finally, note that the number of additional Majorana 
fermions we need does not depend on the number of phys- 
ical fermionic modes at each site: we only need one pair 
of Majorana fermions for all terms in the Hamiltonian 
connecting a particular pair of sites that are non local in 
the qubit representation. 
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Appendix E: Local Unitaries for Discrete Weyl 
Fermions in Three Dimensions 



Appendix F: Infinite Systems 



Here we will see that we can decompose the discrete 
evolution operator U = T1T2T3 into a product of local 
unitaries. Each T, is a conditional shift operator in the 
ith spatial direction that can be rewritten as a product 
of local swap operations. 

It will be convenient to change notation slightly, so 
that 



(E.l) 



Then T 3 , the conditional shift in the z direction, is 

it 



(E.2) 



where e z T — (0,0,1). But we already know how to 
rewrite this in terms of local swap operations, which we 
did in the one dimensional case in equations |IV.11| and 
|IV.12| So T3 is equivalent to applying 



at each ft, followed by 



(E.3) 



(E.4) 



at each ft. Similarly, T% and T2 are conditional shifts de- 
pending on the internal degree of freedom of the particle. 
We define 



so that 



(E.5) 



(E.6) 



where e^ = (1,0,0). And this has a decomposition in 
terms of local swaps also: 

i>(ft,U) ^ V>(fJ-<r*,t„)> ( EJ ) 
at each n, followed by 

%,t.) ^ ( E - 8 ) 
at each n. An analogous decomposition holds for T 2 , with 



(E.9) 



Note that z/>L + * creates a fermion with spin up in the y 
direction. 



Here we will prove our main results when the spatial 
lattice is infinite. Note that, if we are concerned only 
with finite times and finite regions of space, which is the 
case for simulations, we do not need to consider infinite 
lattices, as the evolution is strictly causal. That said, 
the extension to infinite lattices emphasizes the fact that 
causal fermions are analogues of quantum cellular au- 
tomata, and that theorem [l] is a fermionic analogue of 
the main result of [fH [7j , as that result is proved for sys- 
tems on infinite lattices. 

When dealing with quantum systems composed of in- 
finitely many subsystems, it is not clear at first glance 
what we should take as our Hilbert space. This is be- 
cause an infinite tensor product of Hilbert spaces is not 
sensible. One way to get around this is to take local oper- 
ators to be fundamental and represent them by elements 
of an abstract algebra. Then define states as functionals 
of the elements of the algebra. This is the C*-algebra 
approach [2"4] . 

We will now give the precise definition of a C*-algebra. 
First, a complex algebra is a complex vector space with 
a product operation that is associative and distributive 
over addition. 

Definition 2. A C*-algebra A is a complex algebra with 
a norm \\ ■ \\, in which it is complete, and an anti-linear 
map A — > A* , with the following properties: 

1. (AB)* = B* A* 

2. \\AB\\ < \\A\\\\B\\ 



3. \\A* 



\A\\ 



I \\A*A\\ = \\A\\\ 

A simple example of a C*-algebra is A4 n (C), the set of 
n x n complex matrices, where the norm is the spectral 
norm (the largest singular value) and the * operation is 
the hermitian conjugate. Because it is finite dimensional, 
this example misses out on the subtleties associated with 
infinite dimensional vector spaces. 

We assume that all C*-algebras we consider have an 
identity, denoted /. 

Next we define states on the C*-algebra. 

Definition 3. A state on a C*-algebra A is a linear 
functional p that is positive, meaning p(B* B) > VB £ 
A, and normalized, meaning p(I) = 1. 

For an infinite spin chain, a simple example of a state 
is all spins pointing up in the x-direction. In the finite di- 
mensional case, for any state p there is a density operator 
a such that p(A) = tr[aA] for any A in the C*-algebra. 

Next, we require that the evolution is an invertible 
map on the C*-algebra that preserves its structure. Such 
a map is called an automorphism. 

Definition 4. An automorphism a is an invertible linear 
map on the C*-algebra that satisfies 
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1. a{A)a{B) = a(AB) 

2. a(A*) = a(A)* 

3. HA)|| = ||A||. 

The first property implies that the dynamics preserve 
commutation or anticommutation relations. An example 
of an automorphism is A — > U*AU for any U in the C*- 
algebra satisfying UU* = U*U = I. 

As we are working in a discrete spacetime picture, 
there is a subalgebra associated to every spatial point. 
We can define a notion of causality for these systems 
that is a natural extension of the definitions we had in 
the finite case. 

Definition 5. An automorphism a, is causal if, for any 
x and any A localized on x, a(A) is localized on the neigh- 
bourhood of x. 

So, if an automorphism is causal, then observables on 
x cannot spread by more than one step in each direc- 
tion (i.e. the size of the neighbourhood of x) after every 
timestep. 

1. Quantum Lattice Systems and QCA 

The C*-algebra for a quantum lattice system is defined 
by associating elements of the algebra to finite regions of 
the lattice, with the property that elements associated to 
A and A' commute if An A' — 0. Furthermore, the set of 
all elements associated to a finite A is isomorphic to the 
set of operators we get by assigning finite dimensional 
Hilbert spaces to the systems in A. Also, the norm of 
elements in a finite region A is just the operator norm on 
the corresponding operators in the Hilbert space picture. 
So locally the C*-algebra looks like a finite quantum sys- 
tem. For the example of a line of qubits, we have that the 
algebra associated to each site is equivalent to A^C). 

It is a useful result [H] that to specify a state we need 
only specify a family of states p\ on every finite region 
A, with the consistency condition that, if A C A', 

p A (A) = p A ,(A), (F.l) 

where A is an element of the algebra associated to region 
A. 

The precise definition of a Quantum Cellular Automa- 
ton is as follows. 

Definition 6. A Quantum Cellular Automaton is a 
quantum lattice system together with evolution over dis- 
crete timesteps via a causal automorphism. 

2. Fermions 

For fermions, the C*-algebra is generated by objects 
satisfying the canonical anticommutation relations: 

{at, Off} = 5gg ^. 
{as, ay} = 0, 



as in the finite case, but now x e 7L d . We will refer to 
al and creation and annihilation operators even 

though they are no longer operators, rather elements of 
an abstract algebra. Also, we will use a dagger to denote 
the * operation. To simplify notation, we will assume 
that there is only one fermionic mode at each spatial 
point, but all of the following results hold with extra 
degrees of freedom. 

As in the finite case, we say that a fermionic operator is 
localized on a spatial region R if it can be written in terms 
of creation and annihilation operators corresponding only 
to R. 

In section [mi we justified an extra requirement on the 
evolution of fermions: given a system of fermions evolving 
via a unitary U, we can add additional fermionic modes 
whose creation operators anticommute with the origi- 
nal fermion creation operators while commuting with U . 
Here, we make an analogous requirement on the evolu- 
tion. 

Requirement 1. Given a system of fermions evolving 
via a causal automorphism a, we can add additional 
fermionic modes whose creation operators anticommute 
with the original fermion creation operators but are un- 
changed by a. 

Before moving on to the local decomposition for causal 
evolutions, we will prove two useful lemmas. 

Lemma 3. Given a causal automorphism a of fermions 
with annihilation operators as satisfying the above re- 
quirement, a(as) is a linear combination of odd products 
of fermion creation and annihilation operators. 

Proof. We write a(a s ) = A odd + A even , where A odd is the 
sum of all the terms that are products of an odd number 
of creation and annihilation operators and A even is the 
sum of all terms that are products of an even number of 
creation and annihilation operators. 

Our extra requirement above implies that we can add 
a fermionic mode with creation operator b\ which anti- 
commutes with all of the original creation and annihila- 
tion operators while having a (tf) = . But this implies 
that a({b\a s }) = {b\A odd + A even } = 0, which is only 
possible if A even — 0. □ 

Next we give another useful lemma. 

Lemma 4. The inverse of a causal fermionic evolution 
a is also a causal fermionic evolution. 

Proof. Lemmaphells us that a(ag) must be a linear com- 
bination of oddproducts of creation and annihilation op- 
erators. So, since a is causal, 

{a(a$),a s } = 0, 

as rU x n (R3) 

and {a{a^), ag\ = (J 

for all y when x is not in the neighbourhood of y. Since 
or 1 is an automorphism, 

{a^a" 1 ^)} = 0, 

and {al, a (a s )} = 
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for all y when x is not in the neighbourhood of y. But 
x is in the neighbourhood of y implies y is in the neigh- 
bourhood of x. It follows that a~~ 1 (ag) is localized on the 
neighbourhood of x, so a -1 is causal. □ 

We are now ready to decompose causal evolutions into 
products of local unitaries. 



3. Local Decomposition of Causal Evolutions 

We denote the automorphism that swaps a fcrmionic 
mode, with annihilation operator ag, and its copy, 
with annihilation operator bg, by sg. This means that 
Sx(as) = bg and sg(bg) = ag. Note that we denote the 
composition of two automorphisms a\ and oti by a\a2, 
meaning a 1 a 2 (^) = 0:1(0:2 (.A)). 

Theorem 3. Given a system of fermions, evolving via 
a causal evolution a, the evolution of two copies of this 
system via o/3 _1 ; where j3 is equivalent to a but acting on 
the copy system, can be decomposed into local unitaries: 



ap- 1 =]J(s S )l[[ps^- 1 } 



(F.5) 



where sg and f3sfi/3 1 are equivalent to conjugation by 
commuting local unitaries. 



Proof. First, 



where 



n 



(F.7) 



is the global swap. This implies s(3s = a. It follows that 



Furthermore, for fermions, applying j3s$l3~~ l is equiva- 
lent to conjugation by a local fermionic unitary because 
sg(A) = SgASg, where 



Sg = exp[z-(6l - al)(bg - ag)], 



(F.9) 



and so f3sgf3 1 is equivalent to conjugation by 

exp [i^-al)(b'g-ag)], (F.10) 



where b' s — (3(bg), which must be localized within the 
neighbourhood of x because /3 is causal. □ 



4. Representation by Qubits 

We have already seen that Sg is local in the qubit 
representation with a sensible choice of ordering for the 
Jordan- Wigner Transformation. To guarantee that the 
other unitaries are local in the qubit picture, however, we 
may need to apply the trick of adding Majorana fermions 
from appendix Dj We saw that fisgfi~ l is equivalent to 
conjugation by 



exp[z-(&g-4)(?4- 



(F.ll) 



This has the form e~ lHs , where Hg is a self adjoint op- 
erator localized on the neighbourhood of x, so we can 
apply the trick of adding Majorana fermions to get local 
unitaries in the qubit picture. This is always possible be- 
cause b'g = P(bg) must be a finite linear combination of 
odd powers of creation and annihilation operators from 
the neighbourhood of x. 

Note that the number of additional fermions we need 
to add per site is finite for causal evolutions. 

When we considered a finite number of fermion modes, 
we needed the state to be a +1 eigenstate of all the op- 
erators Mgg = im(g^m,(g^ we introduced. Now, in the 
infinite case, we need to construct an analogous state. 
Take any state of the original fermions and extend it to a 
state on the total system of physical and additional Ma- 
jorana fermions p, with the property that all additional 
fermionic modes are empty. Then define the state a by 
a family of states a\ on finite regions A: 



o-a(A) 



with K 



p{K^AK) 
1 , 

feA 



n 



l {y,x) 



(F.12) 



where A is localized on A. Note that there are only 
finitely many rnig g\ for each x because the neighbour- 
hood of x is finite. Also, the ordering of the terms in the 
product is not critical. The state a has the property that 
a(Mrg t g\A) = a(AMrg t g,) = o~{A), (this can be seen from 
equation |G.l | at the start of the following section) so that 



the local fermionic unitaries augmented with the Majo- 
rana fermions are equivalent to the original unitaries on 
this state. Furthermore, the results of a measurement on 
the physical fermions in the state a are the same as those 
from the same measurement on the original state. 

To view this as a QCA, we map the fermionic system 
to a qubit lattice via the Jordan- Wigner Transformation. 
Given any sensible ordering of the infinite lattice sites (for 
example, starting at the origin and spiralling outwards 
filling progressively larger cubes), this is an isomorphism. 
It will map every element of the fermionic C*-algebra to 
an element of the qubit C*-algebra, as the string of Zs 
that arise from the Jordan- Wigner Transformation for 
any ag or m^g^ will be finite. 

So, by adding additional Majorana fermions, we have 
local qubit unitaries Vg implementing the local fermionic 
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unitaries of equation F.ll when in the state a. Note, 



however, that it is not necessarily true that Vg and Vg 
commute when the neighbourhood of x and the neigh- 
bourhood of y overlap. But this does not matter. We can 
implement them in a finite number of steps, where each 
step involves applying Vg unitaries that are localised on 
areas that do not overlap. Note that the order in which 
they are applied does not matter because these opera- 
tors commute when acting on the state a. So we have 
extended theorem [5] to the infinite case. 



Appendix G: Preparing the Majorana state 

We want to prepare a state in the qubit picture that 
is a +1 eigenstate of M/g t g\ = im^^mng^y To do this, 
first notice that, because m fsy) = m \yx) = •"■> 

M (x,y){m(x,y) ~ = ( m (x,y) ~ (G.l) 

It follows that 

II ^|( m (^y) ~ im (y,S))\ n ) ( G - 2 ) 



is a +1 eigenstate of all Majorana pairs M/g^gs, where 
(xy) denotes pairs of sites where we add Majorana 
fermions. The order of the product is irrelevant since 
any order will be a +1 eigenstate of the Mrg t g\ pairs. 
Note that 



II ~m^ m ^,y) ~ im { g S) )\n) 



(Sy) 



(G.3) 



where, m^g) = c\ 3 ^ + c ( s,y), with c^\Q) = 0. This 
state can be created up to a phase by applying the uni- 
taries 



exp|z 



(G.4) 



to \Q). To see this, note that Bfg t g) = 1) which implies 
that 



We want to create the invariant state on qubits, but in 
the qubit representation the creation operators c|- still 
have those awkward strings of Z operators, making this a 
non local unitary. Still, the unitaries in equation |G.4| can 
be implemented efficiently by using a method presented 
in [3] to deal with the strings of Z operators. First, con- 
sider all the qubits that the qubit representation of ct- 
acts on with a Z. We can map the parity of these qubits 
to a flag qubit, which means that a single Z acting on the 
flag qubit has the same effect as the string of Zs applied 
to the other qubits. For example, with € {0, 1}, 



Za...Z n \r ...r n )\22 r j mod 2 ) = 

3=0 
n 

(_l)5:?=o^ m ° d2 | ro ... r „)|^ rj m od 2) = (G.6) 

3=0 

71 

|r — r n )Z\ y^r 3 - mod 2). 

3=0 



So, after preparing flag qubits for c|-,^ and c|-~^ the 



qubit unitary in (G.4 1 is equivalent to a unitary on four 



qubits. After this step, we need to reverse the operation 
preparing the flag qubits, but this and the original flag 
preparation can be done using only n two qubit unitaries, 
where n is the number of qubits we count the parity of. 
For example, if n — 2, we need only two steps: 



\nr 2 )\0) f -> |rir 2 )|n)/ -> 
|rir 2 )|(r 1 + r 2 ) mod 2} f , 



(G.l) 



where g {0, 1} and the subscript / denotes the flag 
qubit. There is a constant number of qubits associated to 
each site because the number of Majorana fermion pairs 
we introduce per site does not grow with Af, so counting 
the parity takes 0(j\f) two qubit operations. Therefore, 
preparing the +1 eigenstate of all Majorana terms can 
be done in time 0(Af 2 ). 

Note that we could also use this method to apply the 
Us^x^Ug gates without introducing additional Majo- 
rana fermions, but this would mean an overhead of 0(j\f) 
for each gate. 

Similarly, to create initial physical fermions states (es- 
sentially, by applying at), we can use the same method 
to deal with strings of Zs appearing in the qubit picture, 
which will add an overhead of 0(j\f) steps. 



